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u,hlS  studyrwe  use  the  National  Centers  for  Environmental  Prediction 
(NCEP)  monthly  sea  surface  temperature  (SST)  fields  (1982-1994)  to  investigate 
the  temporal  and  spatial  variabilities  of  the  South  China  Sea  (SCS)  warm/cool 
anomalies.  Three  steps  of  analysis  were  performed  on  the  data  set:  ensemble 

mean  (T),  composite_analysis  to  obtain  the  monthly  mean  anomaly  relative  to 
the  ensemble  mean  (T),  and  empirical  orthogonal  function  (EOF)  analysis  on 

the  residue  data  relative  to  (T)  +(T).  The  ensemble  mean  SST  field  ( T )  has  a 
rather  weak  horizontal  gradient:  29°C  near  the  Borneo  coast  to  25°-26°C  near  the 
southeast~-China  coast-  Two  areas  of  evident  SST  anomalies  were  found  in  the 
monthly  T  variation:  west  of  Borneo-Palawan  Islands  (WBP)  and  southeast  of  the 
southern  Vietnam  coast  (SVC).  Four  patterns,  monsoon  and  transition  each  with 
two  out-of-phase  structures,  were  found.  During  the  spring-to-summer  transition 
(March  to  May)  the  warm  anomaly  is  formed  in  the  northern  SCS  with  T  >  1  8°C 
located  at  112o-119°30'E,  15°-19°30'N.  During  the  fall-to-winter  transition  (October 
to  November)  the  northern  SCS  (north  of  12°N)  cool  anomaly  is  formed  in  November 
with  T  <  -0.6°C  located  at  108M15°E,  13°-20°N.  We  performed  an  EOF  analysis 
on  the  residue  data  relative  to  T  +  f  in  order  to  obtain  transient  and  interannual 
variations  of  the  SST  fields.  EOF1  accounts  for  47%  of  the  variance  and  represents 
the  northern  SCS  warm/cool  anomaly  pattern.  EOF2  accounts  for  14%  of  the 
variance  and  represents  the  southern  SCS  dipole  pattern.  Strong  northern  SCS 
warm  anomaly  (1°C  warmer)  appears  during  October-November  1987  and  January- 
February  1988  and  strong  northern  SCS  cool  anomaly  (1°C  cooler)  occurs  during 
March  1986  and  November  1992.  Furthermore,  a  strong  cross  correlation  between 
wind  stress  curl  and  SST  anomalies,  computed  from  the  European  Centre  for 
Medium-Range  Weather  Forecast  analyzed  wind  stress  data  and  the  NCEP  SST 
data  for  different  lags,  shows  the  existence  of  an  air-sea  feedback  mechanism  in  the 
bCb  deep  basin. 


1.  Introduction 

The  South  China  Sea  (SCS)  has  a  bottom  topography 
(Figure  1)  that  makes  it  a  unique  semienclosed  ocean 
basin  that  is  overlaid  by  a  pronounced  monsoon  sur¬ 
face  wind.  Extended  continental  shelves  (less  than  200 
m  deep)  are  found  on  the  northern  and  the  southwest¬ 
ern  parts,  while  steep  slopes  with  almost  no  shelves  are 
found  in  the  eastern  part  of  SCS.  The  deepest  water  is 
confined  to  a  bowl-type  trench.  The  maximum  depth  is 
around  4700  m. 
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Based  on  limited  data  sets,  both  cool  and  warm 
anomalies  were  detected  in  SCS.  Dale  [1956]  reported  a 
cool  anomaly  off  the  central  Vietnamese  coast  in  sum¬ 
mer.  Nitani  [1970]  found  a  cool  anomaly  located  at  the 
northwest  of  Luzon.  Reports  from  the  South  China  Sea 
Institute  of  Oceanology  [1985]  indicate  that  in  the  cen¬ 
tral  South  China  Sea,  a  warm  anomaly  appears  in  both 
summer  and  winter  but  closer  to  Vietnam  in  summer 
at  the  surface.  Recently,  a  warm  anomaly  was  reported 
in  the  central  SCS  during  the  late  spring  season  [Chu 
and  Chang,  1995a,  b,  1997;  Chu  et  ai,  1996,  1997]  and 
a  cool  anomaly  was  detected  in  the  central  SCS  during 
December  29,  1993  to  January  5,  1994  from  the  analy¬ 
sis  of  TOPEX/  POSEIDON  data  [5oon5  et  ai,  1995], 
What  are  the  temporal  and  spatial  variabilities  of  sea 
surface  temperature  (SST)  anomaly?  We  used  the  Na¬ 
tional  Centers  for  Environmental  Prediction  (NCEP) 
SST  data  set  to  investigate  these  problems. 
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Figure  1.  Geography  and  isobaths  showing  the  bottom  topography  of  the 
South  China  Sea. 


The  SST  data  (1982-1994)  for  this  study  were  ob¬ 
tained  from  the  NCEP  monthly  real-time  SST  global 
data  set,  which  was  established  by  blending  in  situ 
and  satellite  SST  data  [Reynolds,  1988;  Reynolds  and 
Marsico,  1993].  The  monthly  optimum  interpolation 
(01)  fields  are  derived  from  linear  interpolation  of  the 
weekly  01  fields  to  daily  fields  then  averaging  the  daily 
values  over  a  month  [Reynolds  and  Marsico,  1993]. 
The  monthly  fields  have  the  same  spatial  resolution 
(1°  x  1°)  as  the  weekly  fields.  Before  the  analysis  is 
completed,  the  satellite  data  are  adjusted  for  biases  us¬ 
ing  the  method  of  Reynolds  [1988]  and  Reynolds  and 
Marsico  [1993].  A  description  of  the  OI  analysis  is  given 
by  Reynolds  and  Smith  [1994].  The  bias  correction  im¬ 
proves  the  large-scale  accuracy  of  the  01.  The  data  pro¬ 
cessing  procedure  described  here  was  done  by  NCEP. 

In  situ  SST  data  were  obtained  from  the  Comprehen¬ 
sive  Ocean-Atmosphere  Data  Set  for  the  period  1981- 
1989,  and  from  radio  message  carried  on  the  Global 
Telecommunication  System  for  the  period  1990-present. 


The  satellite  SST  data  were  obtained  from  analysis  of 
National  Environmental  Satellite  Data  and  Information 
Service  data  by  the  University  of  Miami’s  Rosentiel 
School  of  Marine  and  Atmospheric  Sciences.  We  use  the 
monthly  SST  for  the  period  1981-1994  for  this  study. 

2.  Seasonal  Variation  of  the  SCS  Wind 
Stress 

SCS  experiences  two  monsoons,  winter  and  summer, 
every  year.  During  the  winter  monsoon  season,  a  cold 
northeast  wind  blows  over  SCS  (Figure  2a)  as  a  re¬ 
sult  of  the  Siberian  high  pressure  system  located  over 
the  east  Asian  continent.  Radiative  cooling  and  persis¬ 
tent  cold  air  advection  maintain  cold  air  over  SCS.  The 
northeast-southwest  oriented  jet  stream  is  positioned  at 
the  central  SCS.  Such  a  typical  winter  monsoon  pattern 
lasts  nearly  6  months  (November  to  April).  During  the 
summer  monsoon  season,  a  warm  and  weaker  southwest 
wind  blows  over  SCS  (Figure  2b).  Such  a  typical  sum- 
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Figure  2.  Mean  surface  wind  stress  vectors  in  the  South  China  Sea  for  (a) 
December  and  (b)  June  (computed  from  the  ECMWF  data  set). 


mer  monsoon  pattern  lasts  nearly  4  months  (mid-May 
to  mid-September).  The  mean  surface  wind  stress  over 
SCS  is  nearly  0.2  N/m2  and  reaches  0.3  N/m2  in  the 
central  portion  in  December  (Figure  2a)  and  is  nearly 
0.1  N/m2  in  June  (Figure  2b). 

Using  the  European  Centre  for  Medium- Range  Weather 
Forecast  (ECMWF)  analyzed  wind  data,  Trenberth  et 
al.  [1989]  calculated  global  monthly  mean  wind  stress 
curl  C  on  a  2.5°  x  2.5°  grid  for  the  period  1982-1989. 
The  drag  coefficient  was  corrected  for  atmospheric  sta¬ 
bility,  with  sea-air  temperature  difference  and  the  rela¬ 
tive  humidity  taken  to  be  mean  monthly  climatological 
values.  The  ensemble  mean  (1982-1989)  C  field  (Figure 
3a)  shows  a  dipole  pattern:  cyclonic  curls  over  the  east¬ 
ern  SCS  and  anticyclonic  curls  over  the  western  SCS. 
This  coincides  with  the  prevailing  northeast  winds  over 
the  central  SCS. 

Four  patterns,  monsoon  and  transition  each  with  two 
out-of-phase  structures,  were  found  (Figures  3b-3e). 
The  monsoon  pattern  is  featured  by  a  northwest  an¬ 
ticyclonic  (cyclonic)  curl  and  a  southeast  cyclonic  (an¬ 
ticyclonic)  curl  in  winter  (summer).  The  transition  pat¬ 
tern  is  featured  by  a  southern  SCS  cyclonic/anticyclonic 
dipole. 

1.  Fall-to-winter  transition  pattern  shows  a  southern 
SCS  (south  of  12°N)  dipole  curl  structure  located  at  the 
west  of  Borneo-Palawan  Islands  (WBP)  and  the  south¬ 
east  of  the  southern  Vietnam  coast  (SVC).  The  typical 
pattern  is  featured  by  a  WBP  anticyclonic  curl  and  a 
SVC  cyclonic  curl  (Figure  3b).  The  strength  of  the 
curl  is  quite  weak,  and  the  maximum  absolute  value  is 
1.2xlO-7N/m3  for  the  cyclonic  curl  and  0.6x  10-7N/m3 
for  the  anticyclonic  curl. 


2.  Winter  monsoon  pattern:  During  the  winter  mon¬ 
soon  season,  a  northeast-to-southwest  oriented  line  sep¬ 
arates  the  SCS  into  two  parts:  the  western  part  un¬ 
der  anticyclonic  wind  stress  curl  and  the  eastern  part 
is  under  cyclonic  wind  stress  curl,  that  is,  the  AC-C 
pattern  (Figure  3c).  The  maximum  absolute  values  of 
both  curls  reach  4xl0_7N/m3.  Such  a  pattern  persists 
from  November  to  March. 

3.  Spring-to-summer  transition  pattern  shows  an 
opposite  dipole  pattern  to  the  fall-to-winter  transition 
period.  The  dipole  curl  is  characterized  by  a  WBP  cy¬ 
clonic  curl  and  a  SVC  anticyclonic  curl  (Figure  3d). 
The  strength  of  the  dipole  curl  is  much  weaker  than  the 
winter  monsoon  period. 

4.  Summer  monsoon  pattern  is  180°  out-of- phase 
with  the  winter  monsoon  pattern.  A  northeast- to-  south¬ 
west  oriented  line  separates  the  SCS  into  two  parts:  the 
western  part  under  cyclonic  wind  stress  curl  and  the 
eastern  part  is  under  anticyclonic  wind  stress  curl  (Fig¬ 
ure  3e).  The  strength  of  the  wind  stress  curl  during 
summer  monsoon  season  is  about  half  that  during  win¬ 
ter  monsoon  season. 

3.  Composite  Analysis  of  SST 

3.1.  Ensemble  Mean  SST  Field 

We  now  examine  the  data  to  see  if  we  can  obtain  a  warm 
anomaly  signal.  SST  is  represented  by  T(x,-,  yj,  Xk,  ti), 
where  (x,,  yj)  is  the  horizontal  grids,  \k  =  1982,  1983, 

...  1994  is  the  time  sequence  in  years,  and  ti  =1,  2,  ... 

12  is  the  monthly  sequence  within  a  year.  Before  in¬ 
vestigating  the  monthly  variation  of  SST,  we  define  the 
following  two  temporal  averages: 


Latitude  (N)  Latitude  (N) 
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Figure  3.  Surface  wind  stress  curl:  (a)  ensemble  mean,  (b)  September,  (c) 
January,  (d)  April,  and  (e)  July  (computed  from  Trenberth  et  al.  1989).  The 
unit  is  10-8  N/m3. 


1  1994 

T(xi,yj,t,)  =  —  ^  T(xi,yj,Xk,ti), 

*  k= 1982 

A*  =  13  (1994-  1982  +  1).  (1) 

which  is  the  long-term  mean  value  for  the  month  t/,  and 


=  1  12  - 

T(xi>  Vj)  =  TZ  J2  ^  (2) 

£  i=i 

which  is  the  ensemble  mecin.  The  ensemble  mean  (1982- 
1994)  SST  field  over  the  SCS  (Figure  4)  shows  a  pattern 
of  northeast-southwest  oriented  isotherms  with  a  pos- 
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Figure  3.  (Continued) 


itive  temperature  gradient  toward  southeast  near  the 
equator.  The  ensemble  mean  has  a  rather  weak  horizon¬ 
tal  temperature  gradient,  decreasing  from  29°C  near  the 
Borneo  coast  to  25°C  near  the  southeast  China  coast. 

3.2.  Mean  Monthly  SST  Anomalies  (T) 

The  long-term_monthly  mean  values  relative  to  the 
ensemble  mean,  T(xi,yj), 

T(xi,yj,U)  =  T[xi,yj,t,)-T(titVj)  (3) 

leads  to  the  composite  features  of  the  monthly  mean 
SST  anomalies,  which  have  the  following  features  (Fig¬ 
ure  5). 

1.  The  typical  winter  (December  to  February)  T  pat¬ 
tern  contains  northeast-to-southwest  oriented  isotherms 
in  the  northern  SCS  (north  of  12°N),  a  WBP  high 
anomaly  (T  >  -1.2°C  in  January),  and  a  SVC  low 

anomaly  (T  <  -2.4°C  in  January)  in  the  southern  SCS 
(south  of  12°N). 

2.  Spring-to-summer  transition  (March  to  May)  pat¬ 
tern  shows  the  northward  expansion  of  the  WBP  high 
anomaly  and  the  formation  of  a  northern  SCS  (north  of 
12°N)  high  anomaly.  The  isoline  of  T  -  1°C  encloses 
almost  all  the  SCS  in  May.  The  high  anomaly  with 
f  >  1.8°C  is  located  at  112°-119°30'E,  15°-19°30'N.  A 
small  low  anomaly  forms  near  WBP. 

3.  Summer  (June  to  September)  T  pattern  is  char¬ 
acterized  by  northeast-to-southwest  oriented  isotherms 
in  the  northern  SCS  (north  of  12°N),  a  WBP  low 
anomaly  ( T  <  0.4°C  in  July),  and  a  SVC  high  anomaly 


(T  >  1.0°C  in  July)  in  the  southern  SCS.  The  summer 
T  pattern  is  out-of-phase  of  the  winter  pattern. 

4.  Fall-to-winter  transition  (October  to  November) 
pattern  shows  the  northward  expansion  of  a  WBP  low 
anomaly  and  the  formation  of  a  northern  SCS  (north 
of  12°N)  low  anomaly  in  November.  The  low  anomaly 
with  f  <  — 0.6°C  is  located  at  108°-115°E,  13°-20°N 
(out-of-phase  of  the  spring-to-summer  transition  pat¬ 
tern).  A  small  high  anomaly  forms  near  WBP. 

These  features  can  be  easily  seen  from  a  seasonal 
mean  SST  anomaly  fields  (Figure  6):  During  winter 
(December  to  February)  in  the  north  SCS  (north  of 
12°N)  isotherms  are  nearly  parallel  northeast-southwest 
oriented  curves  with  T  decreasing  northward,  and  in  the 
south  SCS  (south  of  12°N)  there  exists  a  WBP  high 
anomaly  and  a  SVC  low  anomaly  (Figure  6a).  Dur¬ 
ing  spring-to-summer  transition,  the  SVC  low  anomaly 
weakens,  and  the  high  anomaly  (original  winter  WBP 
high  anomaly)  moves  toward  west,  expands  northward, 
and  finally  forms  a  large  SCS  high  anomaly.  A  WBP 
low  anomaly  forms  (Figure  6b).  After  T  forms  the  typ¬ 
ical  pattern  shown  in  Figure  6b,  the  SCS  high  anomaly 
starts  to  weaken  and  retreats  toward  SVC.  The  north 
SCS  has  nearly  parallel  isotherms  with  T  increases 
northward,  and  the  WBP  low  anomaly  strengthens. 
Gradually,  it  shifts  to  the  summer  pattern  (Figure  6c), 
which  features  a  WBP  low  anomaly  and  a  SVC  high 
anomaly  in  the  south  SCS.  Such  a  summer  pattern  lasts 
around  4  months  (June  to  September).  At  the  end 
of  summer,  the  low  anomaly  (original  summer  WBP 
low  anomaly)  moves  westward  and  expands  northward, 

ensemble  mean  of  SST  field 
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Figure  4.  The  ensemble  mean  of  the  SST  field  during 
1982-94.  6 
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Figure  5.  Monthly  mean  SST  anomalies  relative  to  the  ensemble  mean  from 
the  NCEP  data  set. 


forming  a  large  SCS  low  anomaly.  A  WBP  warm 
anomaly  forms  (Figure  6d).  After  T  has  the  typical  pat¬ 
tern  shown  in  Figure  6d,  the  SCS  low  anomaly  starts 
to  weaken  and  retreats  to  SVC.  The  north  SCS  has 
nearly  parallel  isotherms,  and  the  WBP  warm  anomaly 
strengthens.  Gradually,  it  shifts  to  the  winter  pattern 
(Figure  6a). 


3.3.  Seasonal  Phase  Shift  Between  (  and  T 

In  order  to  show  the  phase  connection  between  the 
mean  wind  stress  curl  and  the  SST  anomalies,  we  pur¬ 
posely  put  the  September  (fall)  mean  C  as  Figure  3b. 
Comparing  Figure  3  to  Figure  6,  we  found  some  asso¬ 
ciation  between  the  cyclonic  (anticyclonic)  wind  stress 
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Figure  5.  (continued) 
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Figure  5.  (continued) 


4.  Empirical  Orthogonal  Function 
Analysis  of  SST 

The  composite  analysis  shows  the  mean  seasonal  vari¬ 
ation  of  the  SST  anomaly.  Next  we  use  empirical  or¬ 
thogonal  function  (EOF)  analysis  to  investigate  the 
spatial  and  temporal  variabilities  of  the  SCS  high/low 
anomalies. 


4.1.  EOFs  for  the  SST  Anomalies 

The  SST  synoptic  anomalies  obtained  by 

T(xi,yj,Xk,ti)  =T(xi,yj,xk,ti)  -T(xi}yjti)  (4) 

are  rearranged  into  a  NxP  matrix,  T(rn,tp),  n  =  1,  2, 
...,  N;  and  p  =  1,  2,  ...,  P.  Here  P  =  156  is  the  total 
number  of  time  points  used  for  computing  the  covari- 
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Figure  6.  Seasonal  mean  SST  anomaly  (relative  to  ensemble  mean)  pattern 
r  (a)  winter  (December  to  February),  (b)  spring  (March  to  May),  (c)  summer 
une  to  September),  and  (d)  fall  (October  to  November). 
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Figure  7.  The  first  two  EOF  modes  (unit  0.01°C):  (a)  EOF1  and  (b)  EOF2. 
The  solid  (dashed)  lines  indicate  positive  (negative)  values. 


covariance  matrix  is  calculated  by 

flu  R-u  R\n 

„  _  Ri\  R22  •••  R2N 


L  Rni  Rn 2  •••  Rnn  J 

Rnm  =  Jj^f(rn,t,)f(rm,tp)t  N  =  546,  P  =  156 

(5) 

where  n  and  m  (1,  2,  N)  denote  the  grid  locations. 
The  diagonal  elements  of  the  covariance  matrix  Rnn  (n 
=  1,  2,  ...  N)  are  the  variance  at  location  rn.  The  off- 
diagonal  elements  are  the  covariance  with  spatial  lag 
equal  to  the  difference  between  the  row  and  column 
indices.  This  symmetric  matrix  has  N  real  eigenvalues 
\a,  and  eigenvectors  <j>a(rj),  such  that 

N 

Rij<f>a(Tj)  =  \a</>a{ri)y  1  =  1»  2,  ...N  (6) 

j  =  l 

The  eigenvectors  <j>\,  <j>2, ...,  <j>^  are  called  empirical  or¬ 
thogonal  functions.  Each  <j>a  is  a  546-point  (21x26  grid 
in  this  study)  distribution  of  SST  anomaly  pattern.  The 
eigenvalues,  A0(a  =  1, 2, ...,  N),  are  all  positive,  and  the 
summation  of  them,  £  \af  equals  the  total  variance. 
Therefore  A0  is  considered  as  the  portion  of  total  vari¬ 
ance  “explained”  by  the  EOF  <f>a.  It  is  convenient  to 
label  the  eigenfunctions  <t>a  so  that  the  eigenvalues  are 
in  descending  order,  that  is, 


Ai  >  A2  >  A3  >  ...  (7) 

The  data  matrix  T(rn,  tp)  is  thus  approximately  writ¬ 
ten  by 

T(rnJp)  =  Y,PCa(tp)<l>a(Tn)  (8) 

a 

where  PCa{tp)  is  the  principal  component  with  a  unit 
of  °C  and  a  size  of  P,  representing  the  temporal  varia¬ 
tion  of  the  associated  spatial  pattern  described  by  EOF 
4>a(rj). 

4.2.  Principal  EOF  Modes 

Spatial  and  temporal  variabilities  can  be  quantita¬ 
tively  investigated  with  the  method  of  EOF  analysis 
and  a  number  of  its  generalized  forms.  In  order  to  de¬ 
lineate  the  major  modes  of  13  years  (1981-1993),  we 
perform  a  conventional  EOF  analysis  and  obtain  the 
first  six  leading  EOFs,  which  are  able  to  account  for 


Table  1.  Variances  of  the  First  Six  Leading  EOFs 


EOF 

Variance 

Cumulative  Variance 

1 

0.470 

0.470 

2 

0.140 

0.610 

3 

0.093 

0.703 

4 

0.058 

0.761 

5 

0.024 

0.785 

6 

0.017 

0.802 
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Figure  8.  Time  series  of  PCI  for  1982-1994.  Here  two  solid  lines  are  PCI  = 
16.67  and  PCI  =  -16.67,  and  two  dotted  lines  are  PCI  =  8.83  and  PCI  =  -8.83. 


almost  80.1%  of  the  total  variance  during  the  seasons 
(Table  1).  When  combined,  EOF1  and  EOF2  explain 
over  half  (61%)  of  the  total  variance.  Each  EOF  mode 
is  normalized  so  that  its  total  spatial  variance  is  equal 
to  unity.  The  patterns  of  the  first  six  EOFs  are  deemed 
adequate  to  explain  the  SCS  SST  spatial  variabilities 
(Figure  7). 

The  EOF1  mode  (Figure  7a)  shows  a  large  warm/cool 
anomaly  feature  in  the  northern  SCS  (north  of  12°N), 
which  may  reflect  some  nonseasonal  process.  The  neg¬ 
ative  values  of  EOF1  indicate  a  northern  SCS  warm 
anomaly  if  PCi  <  0  and  a  northern  SCS  cool  anomaly 
if  PC\  >  0.  Thus,  we  call  EOF1  the  northern  SCS 
warm/cool  anomaly  mode. 

The  EOF2  mode  (Figure  7b)  shows  a  strong  north-to- 
south  contrast  across  the  SCS,  which  probably  reflects 
the  interannual  fluctuation  of  monsoons.  If  PCi  is  pos¬ 


itive,  the  EOF2  mode  shows  a  decrease  of  SST  with 
latitude.  If  PCi  is  negative,  the  EOF2  mode  shows 
an  increase  of  SST  with  latitude.  We  call  EOF2  the 
north-south  contrast  mode. 

4.3.  Temporal  Variabilities 

Northern  SCS  warm/cool  anomaly  mode.  The 
first  principal  component  PC1(tp )  for  1982-1994  shows 
the  temporal  variability  of  the  northern  SCS  warm/cool 
anomaly  (Figure  8).  Since  the  EOF1  mode  ^(rn) 
is  always  negative  throughout  the  whole  SCS  (Fig¬ 
ure  7a),  positive  values  of  PCi(tp)  correspond  to  neg¬ 
ative  SST  anomalies  (north  SCS  cool  anomaly)  and 
negative  values  of  PCi(tp)  correspond  to  positive  SST 
anomalies  (north  SCS  warm  anomaly).  We  see  in¬ 
terannual  variabilities  in  PCU  over  periods  of  2  to  5 
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Figure  8.  (continued) 


years.  The  maximum  value  of  PC\  is  20,  appearing 
in  November  1992,  which  indicates  a  strong  northern 
SCS  cool  anomaly  with  mean  monthly  SST  anomaly  - 
1.2°C  [—20  x  (— 0.06°C)]  in  that  month.  The  minimum 
value  of  PC\  is  -22,  appearing  in  February  1988,  which 
indicates  a  strong  northern  SCS  warm  anomaly  with 
mean  monthly  SST  anomaly  1.3°C  [—22  x  (— 0.06°C)] 
in  that  month.  If  the  isoline  of  (-0.06°C)  is  treated  as 
the  boundary  of  the  northern  SCS  warm/cool  anomaly 
(Figure  7a),  and  if  we  are  interested  in  evident  warm/cool 
anomaly  occurrence  with  SST  anomaly  exceeding  1°C, 
we  found  the  following  criterions:  When  PCi(tp)  > 
16.67,  the  northern  SCS  has  a  cool  anomaly  with  tran¬ 
sient  SST  anomaly  of  -1°C;  When  PCi(tp)  <  —16.67, 
the  northern  SCS  has  a  warm  anomaly  with  transient 
SST  anomaly  of  1°C.  If  we  consider  a  warm  (or  cool) 


anomaly  with  transient  SST  anomaly  of  0.5°C  (or  - 
0.5°C),  the  criterion  becomes  PC\(tp)  >  8.83  (or  PCi(tp) 
<  -8.83).  Four  lines,  PCX  -  ±16.67  (solid)  and  PCi  = 
±8.83  (dashed),  were  drawn  in  Figure  8  for  identifying 
strong  northern  SCS  transient  SST  anomalies.  Table 
2  lists  the  periods  of  strong  (|AT|  >  1°C)  and  evident 
(|AT|  >  0.5°C)  transient  SST  anomalies.  We  found  two 
strong  warm  anomaly  (AT  >  1°C)  periods  (October- 
November  1987  and  January- February  1988)  and  two 
strong  cool  anomaly  (AT  <  -1°C)  periods  (March  1986 
and  November  1992). 

North-south  contrast  mode.  The  second  princi¬ 
pal  component  PCitf)  for  1982-1994  is  shown  in  Fig¬ 
ure  9.  The  EOF2  mode  (Figure  7b)  02 (r„ ),  reveals 
the  north-south  contrast  structure.  Identification  of  the 
pattern  (cool  northern  SCS  and  warm  southern  SCS  or 
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Table  2.  Northern  SCS  Warm  and  Cool  Anomaly  Periods 


Anomaly  Period 

Above  1°  C 

Above  0.5°  C 

Warm 

Oct. -Nov.  1987,  Jan. -Feb.  1988 

Jan.-Dee.  1982,  Aug.-Oct.  1983 
Jul.  1987 
Oct.  1987-March  1988 
May  1988,  July-Sept.  1988 
June-July  1993,  Feb.  1994 
April-May  1994,  Dec.  1994 

Cool 

March  1986,  Nov.  1992 

Dec.  1983-Feb.  1984.  Aug.  1984 
July  1985,  Sept.-Oct.  1985 
Jan.-March  1986 
Nov.-Dee.  1988 
Nov.  1991 
Oct.-Nov.  1992 
Feb.  1993 
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Figure  9.  The  same  as  Figure  8  except  for  PC2. 
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Figure  9.  (continued) 


vice  versa)  depends  on  the  sign  of  PCi(t).  PCaft)  >  0 
corresponds  to  the  decrease  of  SST  with  latitude  pat¬ 
tern,  and  PCi(i)  <  0  corresponds  to  the  increase  of  SST 
with  latitude  pattern.  The  maximum  absolute  value  of 
PCi  is  13,  appearing  in  February  1983,  which  indicates 
a  near  2°C  south- to-north  decrease  of  SST.  The  mini¬ 
mum  value  of  PCi  is  -11,  appearing  in  February  1987, 
which  indicates  a  near  1.7°C  north-to-south  decrease  of 
SST. 

5.  An  Air-Sea  Feedback  Scenario 

5.1.  Cross  Correlation  Between  and  T 

We  interpolated  the  wind  stress  curl  (£)  data(2.5°  x  2.5°) 
into  the  SST  grid  (1°  x  1°)  and  chose  the  period  1982- 
1989  for  the  cross  correlation  analysis.  Cross  correlation 
between  £  and  T  is  computed  at  each  grid  point  by 


(<at  -  1) 
(*N  ~  |r|) 


£(  {[c(*.y,0  -C(*>y) 

[r{x,y,t  +  r)  -  T[x,y) 

} 

[c(*,y,<)  -C(z.y)]  \jllt 

T(x,y,t)  -  T(x,  y) 

r 

(9) 

Here  t n  =  96,  the  total  months  of  the  data  sets.  When 
t  >  0,  (9)  represents  correlation  between  r  months’ 
leading  of  wind  stress  curl  anomaly  to  SST  anomaly, 
and  r  <  0  denotes  correlation  between  |r|  months’  lead¬ 
ing  of  SST  anomaly  tc  wind  stress  curl  anomaly. 


5.2.  Two  Types  of  Cross-Correlation 
Coefficients 

The  no-lag  cross  correlation  coefficient  (CCC)  be¬ 
tween  C  and  T  (Figure  10)  shows  (1)  negative  CCC  in 
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Figure  10.  No  lag  cross-correlation  coefficient  between 
the  surface  wind  stress  curl  and  SST  anomalies. 


the  deep  basin  of  SCS  (-0.6,  minimum)  and  (2)  positive 
CCC  in  the  shelf  regions  of  SCS  (0.6,  maximum).  This 
basin  negative  -  shelf  positive  (Bn-Sp)  pattern  continues 
as  T  lags  C  (oceanic  response  to  the  atmosphere)  from 
1  to  3  months  (Figures  lla-llc).  Both  positive  and 


negative  R(^T  strengthens  when  the  lag  begins. 
has  a  minimum  value  of  -0.7  and  a  maximum  value  of 
0.8  as  T  lags  C  for  1  month  (Figure  11a).  As  r  further 
increases  (r  >  2  months),  both  positive  and  negative 
R[t)t  weakens  (Figures  11c  and  lid).  When  T  lags  C 

for  5  or  more  months  (Figures  Ilf  and  llg),  rQ  tends 
to  have  an  out-of-phase  pattern:  basin  positive  -  shelf 
negative  (Bp-Sn)  pattern. 

As  C  lags  T  (atmospheric  response  to  the  ocean)  for 
a  month,  the  Bn-Sp  pattern  continues,  but  the  CCC 
values  are  greatly  reduced  (Figure  12a).  As  r  =  -2 
month,  the  Bn-Sp  pattern  disappears  (Figure  12b).  As 
7"  =  -3  months  (Figure  12c),  the  Bp-Sn  pattern  has 
been  established  and  sustained  for  4  months  (Figures 
12d-12f). 

The  dependence  of  CCC  on  r  (Figure  13)  shows  an 
oscillatory  mode  with  a  strong  annual  signal  at  118°E, 
18°N  (in  the  SCS  deep  basin)  and  at  110°E,  14°N  (in 
the  Vietnam  coast).  For  the  deep  basin  point,  CCC 
has  a  value  of  -0.6  at  r  =  0  (no  time  lag),  decreases 
with  r  to  a  minimum  of  -0.8  at  r  =  1  month  (strongest 
negative  correlation),  and  then  increases  with  t  to  a 
value  near  0  at  r  =  4  months  (no  correlation).  CCC 
increases  from  0  to  0.8  (maximum)  as  the  lag  t  increases 
from  4  to  7  months  and  decreases  again  after  r  =  7 
months.  As  r  decreases  from  0  to  -2  (£  lags  T  from  0  to  2 
months),  CCC  increases  from  -0.6  to  0  (no  correlation). 
As  r  decreases  from  -2  to  -5  (C  lags  T  from  2  to  5 
months),  CCC  increases  from  0  to  0.8  (strong  positive 
correlation).  The  out-of-phase  feature  between  the  two 


Figure  11.  Cross-correlation  coefficient  between  the  surface  wind  stress  curl 
and  SST  anomalies  for  various  SST  lags:  (a)  1  month,  (b)  2  months,  (c)  3 
months,  (d)  4  months,  (e)  5  months,  and  (f)  6  months. 
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Figure  11.  (continued) 


regions  (deep  basin  and  shelf)  indicates  different  physics 
to  be  involved  in  the  air-sea  feedback  processes. 

5.3.  Possible  Air-Sea  Feedback  Mechanism  in 
the  Northern  SCS  Basin 

We  may  suggest  an  air-sea  feedback  mechanism  in 
the  SCS  deep  basin  from  the  evident  cross  correlation 
between  C  and  T  for  varying  lag  r.  The  negative  values 
of  (0<  r  <  3)  in  the  basin  indicate  the  Ekman- 


pumping  effect:  cyclonic  wind  stress  curl  induces  up- 
welling  which  sucks  the  deep  cold  water  to  the  surface, 
and  anticyclonic  wind  stress  curl  generates  downwelling 
which  prevents  the  deep  cold  water  to  the  surface.  The 
cold  (warm)  SST  anomaly  tends  to  produce  high  (low) 
surface  pressure  perturbation  which  leads  to  the  genera¬ 
tion  of  anticyclonic  (cyclonic)  wind  stress  curl  anomaly. 
Thus  the  values  of  #{/•}.  become  positive  in  the  basin 
as  -2  >  r  >  —7.  The  wind  effect  to  generate  the  cen¬ 
tral  SCS  warm/cool  SST  anomalies  has  been  verified 
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Figure  12.  Cross-correlation  coefficient  between  the  surface  wind  stress  curl 
and  SST  anomalies  for  various  wind  stress  curl  lags:  (a)  1  month,  (b)  2  months, 
(c)  3  months,  (d)  4  months,  (e)  5  months,  and  (f)  6  months. 


by  a  numerical  model  (Chu  et  al.,  Generation  of  the 
South  China  Sea  deep  basin  warm  and  cool  anomalies 
by  wind-topography  forcing,  submitted  to  Journal  of 
Geophysical  Research,  1997). 

6.  Conclusions 

The  goal  of  our  study  was  to  detect  whether  warm/cool 
anomaly  exists  in  the  SCS  using  the  National  Meteoro¬ 
logical  Center  (NCEP)  monthly  sea  surface  temperature 


(SST)  fields  (1982-1994).  Composite  and  EOF  analy¬ 
ses  were  used  for  the  study.  We  obtained  the  following 
results  from  this  observational  study: 

1.  The  ensemble  mean  SST  field  (f)  was  established 
with  a  rather  weak  horizontal  gradient  (29° C  near  the 
Borneo  coast  to  26°C  near  the  southeast  China  coast). 
This  result  agreed  quite  well  with  the  early  study  using 
the  Navy’s  Master  Oceanographic  Observation  Data  Set 
[Chu  et  al.,  1996,  1997]. 
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Figure  12.  (continued) 


2.  The  composite  analysis  indicates  features  of  the 
monthly  averaged  SST  anomaly  fields  relative  to  the 
ensemble  mean.  Four  patterns,  monsoon  and  transi¬ 
tion  each  with  two  out-of-phase  structures,  were  found. 
The  monsoon  pattern  features  northeast-to-southwest 
oriented  isotherms  in  the  northern  SCS  and  a  dipole 


structure  in  the  southern  SCS.  This  thermal  dipole  is 
out-of-phase  from  winter  to  summer:  The  dipole  with 
the  WBP  warm  anomaly/SVC  cool  anomaly  is  found 
in  winter  and  with  the  WBP  cool  anomaly/SVC  warm 
anomaly  in  summer.  The  southern  SCS  dipole  fea¬ 
tures  a  WBP  warm  anomaly  (T  >  — 1.2°C  in  January) 


Figure  13.  Cross-correlation  coefficient  between  the  surface  wind  stress  curl 
and  SST  anomalies  for  various  lags  at  118°E,  18°N  (circles,  in  the  central  SCS 
deep  basin)  and  110°E,  14°N  (asterisks,  in  the  Vietnam  shelf). 


20,955 


CHU  ET  AL.:  SOUTH  CHINA  SEA  SURFACE  TEMPERATURE  VARIABILITIES 


and  a  SVC  cool  anomaly  (T  <  — 2.4°C  in  January) 
in  winter  (December  to  February);  and  a  WBP  cool 
anomaly  (T  <  0.4°C  in  July)  and  a  SVC  warm  anomaly 
(T  >  1.0°C  in  July)  in  the  summer  (June  to  Septem¬ 
ber). 

The  transition  pattern  is  characterized  by  the  west¬ 
ward  expansion  of  the  WBP  warm  (cool)  anomaly  and 
the  formation  of  the  SCS  warm  (cool)  anomaly  in  the 
spring-to-summer  (fall-to-winter)  transition.  During 
the  spring-to-summer  transition  (March  to  May),  the 
warm  anomaly  is  formed  in  the  northern  SCS  with 
f  >  1.8°C  located  at  112o-119°30'E,  15°-19°30,N.  Dur¬ 
ing  the  fall-to-winter  transition  (October  to  Novem¬ 
ber),  the  northern  SCS  (north  of  12°N)  cool  anomaly 
is  formed  in  November  with  T  <  -0.6°C  located  at 
108°-115°E,  13°-20°N. 

3.  Four  patterns,  monsoon  and  transition  each  with 
two  out-of-phase  structures,  were  found  in  the  surface 
wind  stress  curl.  The  monsoon  features  a  northwest  an- 
ticyclonic  (cyclonic)  pattern  and  a  southeast  cyclonic 
(anticyclonic)  pattern  in  winter  (summer).  The  transi¬ 
tion  pattern  features  a  southern  SCS  cyclonic/anticyclonic 
dipole.  We  found  associations  between  the  cyclonic 
(anticyclonic)  wind  stress  curl  and  the  warm  (cool) 
anomaly  with  a  seasonal  phase  shift.  For  the  south¬ 
ern  SCS,  the  fall-to-winter  (spring-to-summer)  (  dipole 
pattern  corresponds  well  with  the  winter  (summer)  T 
dipole  pattern.  For  the  central  SCS,  the  winter  (sum¬ 
mer)  £  pattern  corresponds  well  with  the  spring  (fall) 

T  pattern.  We  may  claim  that  the  winter  (summer)  £ 
pattern  causes  the  northward  expansion  of  SVC  warm 
(cool)  anomaly  in  the  spring  (fall). 

4.  EOF  analysis  was  applied  to  the  residual  SST  data 
(relative  to  T  +  T)  to  investigate  the  transient  and  in¬ 
terannual  SST  variabilities.  EOF1  accounts  for  47%  of 
the  variance  and  shows  the  northern  SCS  warm /cool 
anomaly  pattern.  The  strongest  northern  SCS  cool 
anomaly  appears  in  November  1992  with  mean  monthly 
SST  anomaly  -1.2°C,  and  the  strongest  northern  SCS 
warm  anomaly  shows  up  with  mean  monthly  SST  anomaly 
1.3°C  in  November  1987  and  February  1988.  We  also 
found  two  strong  warm  anomaly  (AT  >  1°C)  periods 
(October-November  1987  and  January-February  1988), 
and  two  strong  cool  anomaly  (AT  <  -1°C)  periods 
(March  1986  and  November  1992). 

5.  An  air-sea  feedback  mechanism  was  proposed  from 
the  cross  correlation  analysis  between  £  and  T.  We  will 
use  a  coupled  air-sea  model  to  verify  this  hypothesis. 
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